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Abstract 

The spreading characteristics of jets from 
several asymmetric nozzles are studied in comparison 
to those of an axisymmetric jet, over the Mach number 
(Mj) range of 0.3 to 1.06. The effect of tabs in two 
cases, the axisymmetric nozzle fitted with four tabs 
and a rectangular nozzle fitted with two large tabs, is 
also included in the comparison. Compared to the 
axisymmetric jet, the asymmetric jets spread only 
slightly faster at subsonic conditions, while at super- 
sonic conditions, when screech occurs, they spread 
much faster. Screech profoundly increases the spread- 
ing of all jets. The effect varies in the different stages 
of screech, and the corresponding unsteady flowfield 
characteristics are documented via phase-averaged 
measurement of the fluctuating total pressure. An or- 
ganization and intensification of the azimuthal vortical 
structures under the screeching condition is believed to 
be responsible for the increased spreading. Curiously, 
the jet from a ’lobed mixer’ nozzle spreads much less 
at supersonic conditions compared to all other cases. 
This is due to the absence of screech with this nozzle. 
Jet spreading for the two tab configurations, on the 
other hand, is significantly more than any of the no- 
tab cases. This is true in the subsonic regime, as well 
as in the supersonic regime in spite of the fact that 
screech is essentially eliminated by the tabs. The 
dynamics of the streamwise vortex pairs produced by 
the tabs cause the most efficient jet spreading thus far 
observed in the study. 

Nomenclature 

C f Measured thrust normalized by ideal thrust. Ideal 

thrust is obtained from actual mass flux through 
nozzle, and assuming uniform, convergent flow. 
D Equivalent diameter based on nozzle exit area. 

f s Screech frequency. 


m Mass flux obtained from Pitot probe survey at a 
given x. 

Tn c Mass flux at nozzle exit obtained from P T /p A 
assuming uniform flow. 

M Mach number. 

Mmax Maximum Mach number at a given x. 

M, Jet Mach number assuming fully expanded flow 
for a given P T /p A . 
p A Ambient pressure. 

P T Plenum chamber pressure. 

St Strouhal number of screech, f s D/U E . 

T Jet thrust. 

U Mean streamwise velocity. 

U E Jet velocity at nozzle exit. 

Umax Maximum U at a given x. 

x,y,z Streamwise and transverse coordinates. 

1. Introduction 

Jets from asymmetric nozzles have been known 
to spread and mix faster than their axisymmetric coun- 
terparts. For example. Refs. 1 and 2 had shown that the 
overall spreading of small aspect ratio elliptic jets was 
clearly more than that of an axisymmetric jet. There were 
many studies on jets from rectangular nozzles, 3 " 7 some of 
which also indicated more vigorous entrainment charac- 
teristics. While most of the studies cited above were for 
low Mach number subsonic jets, a similar inference was 
also made for supersonic elliptic jets. 8 The possible impact 
of faster jet spreading on noise reduction was addressed 
in Ref. 9. 

Based on the available data, however, a clear and 
direct comparison of the spreading characteristics of jets 
from the various nozzles has been difficult. Part of the 
difficulty is due to the fact that mixing and spreading of 
a jet can be defined in a number of ways, and different 
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studies have used different parameters for that purpose. 
Many parameters, such as streamwise variation of the jet 
half-velocity-width or the centerline velocity, would only 
provide a partial description and could even be mis- 
leading. Another difficulty stems from differences in the 
operating conditions in the previous studies, especially in 
Mach number. It is well known that mixing layer spread- 
ing decreases with increasing compressiblity at higher 
Mach numbers. 10 Therefore, a comparison of jet spreading 
with different nozzles should be meaningful only when 
done at the same Mach number and preferably based on 
the measurement of the whole flow field. The objective of 
the present study is to conduct such a comparative evalua- 
tion. 

In the continuing effort to increase jet mixing and 
spreading, various other methods, e.g., application of 
vortex generators in the form of ’tabs’, are also under 
investigation. 11 ' 16 A tab is a small protrusion placed at the 
jet nozzle exit which produces counter-rotating streamwise 
vortex pairs. The streamwise vortex pairs, depending on 
their number, sense and strength, can have a profound 
inpact on the jet spreading. The effect will be summa- 
rized in the appropriate context. In the present com- 
parative evaluation, two tab configurations are also 
included. Both configurations, to be described further 
shortly, exhibited significant increase in the jet spread- 
ing. 12 - 14 

The bulk of the data, presented in the following, 
pertain to a fixed subsonic and a fixed supersonic con- 
dition. In addition, jet spreading characteristics are 
evaluated as a function of Mach number for limited cases. 
In particular, the various stages of screech at the super- 
sonic condition, occurring with varying Mach number, 
affect the jet spreading significantly. This is explored for 
the different nozzle geometries. Furthermore, since the 
use of the tabs is accompanied by a thrust penalty which 
is of utmost importance in many applications, 17 ’ 18 thrust 
loss is also measured and analyzed for all the cases. 

2. Experimental Method 

Most of the experiments were carried out in a 
small jet facility. 12 Compressed air supplied through one 
end of a plenum chamber, fitted with flow conditioning 
units, exited through a 3.81 cm diameter ’nipple’ on the 
opposite end. The supply air had a total temperature 
approximately equal to the ambient temperature. The dif- 
ferent nozzles could be attached to the nipple. All asym- 


metric nozzles were machined from solid cylindrical 
blocks of aluminum by electrical discharge machining 
(EDM). The interior along the major and the minor axes 
for each nozzle was contoured according to fourth order 
polynomial fits. The rest of the interior was faired 
smoothly. In all cases the flow always converged, and 
entered and exited the nozzle approximately axially. The 
equivalent diameter (D) was the same for all the asym- 
metric nozzles and was 1.47 cm. The circular nozzle was 
fabricated earlier, and had an exit diameter of 1.27 cm. 
The general interior shape of each nozzle is shown 
schematically in Fig. 1(a). All had ’aid walls’ to facilitate 
easy installation of the tabs at desired locations. 

The asymmetric nozzles included a 3:1 elliptic, 
a 3:1 rectangular and a generic 6-lobed case. The two tab 
configurations were: four equally spaced tabs with the 
circular nozzle and two large tabs spanning the narrow 
edges of the 3:1 rectangular nozzle. Thus, altogether six 
nozzle geometries were considered, and these are shown 
schematically in Fig. 1(b). The tabs used were the ’delta 
tabs’, 12 having triangular shapes with the base on the 
nozzle wall and the apex leaning downstream at 45°. The 



Fig. 1 Schematics of (a) nozzle interior shape and (b) 
exit geometries of four nozzles with and without tabs. 


2 



tab sizes were such that the area blockage due to each was 
about 2% of the nozzle exit area for the circular case 
(total blockage being about 8%), and about 6% for die 
rectangular case (total blockage being about 12%). These 
tab configurations were selected after conducting para- 
metric studies designed to yield a large increase in the 
downstream jet spreading. 12,14 

The mass flux measurements woe made at 
various downstream locations by standard total pressure 
surveys with a 0.76 mm (o.d.) Pitot tube. The measure- 
ments were conducted for a fixed subsonic (M,=0.3) and 
a fixed supersonic (M,= 1.63) condition. Since the nozzles 
were convergent, underexpanded conditions existed in the 
supersonic regime. The jet Mach number M, in that 
regime was defined simply as that which would be 
achieved had the flow expanded fully for a given P T /p A . 
With the supersonic jets, the measurements were done far 
Miftngh downstream so that the flow was subsonic every- 
where. This was done to avoid measurement mors and 
probe interference that are typically encountered in super- 
sonic flows. 

Throughout the data acquisition, the plenum 
pressure ambient pressure were monitored and data 
normalization was done according to current conditions. 

averaging time was allowed to ensure good data 
repeatability. The mass flux results were obtained by in- 
tegration of the Pitot probe survey data over the cross 
sectional plane. It should be noted that there is subjec- 
tivity in the calculation of mass flux. Crow and Cham- 
pagne 19 eloquently stated, "The notions of volume flux 
and entrainment are creatures of theory, in the case of a 
jet, rather than experiment." This is because the integrand 
does not fall off to zero in the surrounding potential flow 
region «nd, "...the volume flux in the induced potential 
flow is infinite." In all previous experiments, therefore, 
various criteria were followed to discriminate the vortical 
flow from the surrounding potential flow. Crow and 
Champagne did this for their hot wire data by plotting the 
integrand, radius times axial velocity, and visually identi- 
fying the potential tails of the flux profiles. Note further- 
more that both hot wire and Pitot probe measurements 
involve errors on the outer periphery of the jet due to 
large turbulence and flow angularity. (With a single hot 
wire, the measured u-profile in the potential tail involves 
slowly diminishing but erroneous positive values pri- 
marily due to radial ent rainme nt velocity. In data obtained 
with a Pitot probe, aligned with the jet axis, small sub- 


ambient values are erroneously read in the tail region, 
again, due to radial entrainment velocity.) 

In view of these difficulties, the criterion fol- 
lowed here simply involved truncating the integration 
where the measured Mach number dropped below 1% of 
die local centerline Mach number. In any case, because of 
the subjectivity noted in the foregoing, it is felt that speci- 
fying a ma rgin of error in the mass flux data is not 
meaningful However, with the given truncation criterion, 
care was falren such that die mass flux data were repeat- 
able within ±4% for the supersonic conditions, the 
repeatability b eing better at upstream locations where the 
velocities were higher. At die subsonic conditions, 
because of smaller dynamic pressures, the data scatter was 
larger and about ±8%. The error in the thrust measure- 
ment was estimated to be well within ±1 %. 

t imiteH data on the vorticity field for an in- 
compressible case are also presented. The measurement 
technique will be briefly described with the data. Phase 
averaged unsteady flowfield data were also obtained for 
die screeching conditions. The periodic pressure fluc- 
tuations were measured as a function of phase over a cross 
seet imml plane of the jet. These measurments were made 
with a 2.34 mm diameter pressure transducer (ENDEVCO 
8507C-50). The probe measured fluctuating total pressure 
with suffici ent frequency response (manufacturer specific- 
ation about 200 kHz). However, the phase-averaged 
measurements could only be done at locations rather close 
to die nr>**le exit where die core flow was still super- 
sonic. Limited surveys showed that the periodicity in the 
flow was essentially lost farther downstream. Thus, these 
measurements involved some probe interference and die 
results should be considered qualitative. However, the 
results quite adequately reveal die mode shapes associated 
with the screech stages. Further details of data acquisition 
and errors will be given with the results. 

3. Results and discussion 

Figure 2 shows an iso-surface of mean velocity 
at the subsonic condition (Mj = 0.3) for die six cases. 
The amplitude represented by die iso-surface is 30% of 
the local maximum velocity; the plots are based on data 
at five x-stations. The two cross-sectional cuts (’surfers’) 
are located at x/D = 2 and 19. A visual impression of the 
jet «Fp«mding for each case can be obtained from these 
data The major axes of the elliptic and die rectangular 
nozzles were aligned with the z-axis. For the lobed 
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Fig. 2 For Mj = 0.3, iso-surfaces of mean velocity, U/U MAX 
range of 2 to 19. 

nozzle, the middle two lobes were aligned with the z-axis. 

For all cases without the tabs the jet cross section be- 
comes round by the end of the measurement range, and no 
significant difference in the spreading can be discerned. 

The two tab cases, in (e) and (f), however, involve 
noticeable increase in the spreading. Note that for the jet 
in (f), as it can be seen from the ’surfer’ on the left, there 
has already been an ’axis switching’ by x/D = 2. For an 
asymmetric jet, the axis switching phenomenon is another 
facet of the spreading characteristics. The dynamics of the 
large scale vortical structures govern both axis switching 
and spreading, as to be discussed further in §3.2. 

Corresponding data for the supersonic condition 
(Mj = 1.63) are shown in Fig. 3. In order to save space, 
the coordinate origin in these data sets have been shown 
at x/D = 8. The ’surfer’ on the right in all cases is at 
x/D = 30, the left ’surfer’ locations are indicated in the 
figure caption. For some cases the measurements were not 
made at x/D = 10 because the core flow was supersonic. 
Clearly, the elliptic, the rectangular, and the tab cases 



= 0.3, for the six cases of Fig. 1(b). Data cover x/D 

result in more spreading. One also finds that unlike at the 
subsonic condition (Fig. 2), the rectangular and the 
elliptic jets (Fig. 3b and 3c) have undergone clear axis 
switching at the supersonic condition. Surprisingly, the 
lobed nozzle involves minimal spreading. These points are 
discussed further in the following. Note that the jet 
spreading in Figs. 3(b), (c) and (f) may appear somewhat 
exaggerated due to the viewing angle. 

The streamwise variation of the maximum 
velocity for the six subsonic cases of Fig. 2 is shown in 
Fig. 4. Consistent with the overall jet spreading seen in 
Fig. 2, one finds the velocity decay to be the fastest for 
the rectangular jet with two tabs. As expected, the slowest 
decay occurs for the circular jet without tab. It is note- 
worthy that the circular jet with four tabs decay faster 
than any of the three asymmetric cases. Similarly, the 
variations of the maximum Mach number for the super- 
sonic cases of Fig. 3 are shown in Fig. 5. Again, consis- 
tent with jet spreading seen in Fig. 3, the Mach number 
decay is found to be the slowest for the lobed nozzle and 
the fastest for the rectangular nozzle with two tabs. 
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Fig. 3 For Mj = 1.63, iso-surfaces of Mach number, M/M max = 0.3, for the six cases. Data cover x/D range of 14-30 
in (a), (c) and (d), and 10-30 in the rest of the cases. Co-ordinate origin is at x/D — 8. 



Fig. 4 For Mj — 0.3, streamwise variations of maxi- 
mum mean velocity for the six cases. 



Kg. 5 For M, = 1.63, streamwise variations of maxi- 
mum Mach number for the six cases. 
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Kg. 6 For Mj = 0.3, streamwise variations of normal- 
ized mass flux for the six cases. 


Normalized mass flux variations with streamwise 
distance for the subsonic cases are shown in Fig. 6. These 
variations correspond to the Mach number decay shown 
in Fig. 4, a larger increase in the former is accompanied 
by a faster decay in the latter. The fluxes for the circular 
jet are the smallest at all x. Corresponding fluxes for the 
elliptic, rectangular and the lobed cases are somewhat 
higher. However, they are not as high as reported in some 
earlier studies. 1 * 2 In ref. 1, for example, the mass flux 
measured at x/D = 5 for a 2:1 elliptic jet was about 55% 
higher than that for a circular jet. The corresponding 
increases for the elliptic and the rectangular jets in Fig. 6 
are less than 10%. A possible reason for this difference 
may trace to differences in initial conditions, and this will 
be discussed further in §3.3. 

The cases with the tabs (Fig. 6) clearly exhibit 
much larger jet spreading. It can be seal from Fig. 2 that 
the cross sections of all the jets, except the one in (f), 
have become essentially round by the last measurement 
station. These have apparently reached the asymptotic, 
self-similar state and, thus, the entrainment rates, given 
by the slopes of the flux curves, are comparable. The 
rate, ((D/ri%).(9m/3x)), based on the farthest two data 
points, turns out to be about 0.27. This is in good 
agreement with previous round jet data. 19 In contrast, the 



fig. 7 For Mj = 1.63, streamwise variations of nor- 
malized mass flux for the six cases. 

rectangular jet with tabs in (f) is still undergoing entrain- 
ment at a higher rate yielding the higher fluxes. One 
might expect that eventually this jet would also attain a 
round cross section. When that happens farther down- 
stream, the entrainment should reach the same asymptotic 
rate as exhibited by the others. (In obtaining the data for 
the tab cases, the area blockage was taken into account. 
That is, the initial flux was calculated by assuming the 
nozzle exit area to be less by the amount of the geometric 
blockage imparted by the tabs. Note that even if this were 
not done, the normalized fluxes would be somewhat lower 
but the trends discussed would still be valid). 

The corresponding mass flux variations for the 
supersonic condition are shown in Fig. 7. A few obser- 
vations can be made, and these will be the basis for 
discussion in the following subsections. First , one finds 
that the circular jet has attained a constant slope, the 
entrainme nt rate being about 0.18. This is much lower 
than that observed in the subsonic case. Second, the 
circular jet with four tabs exhibits a much higher flux in 
the upstream region, indicating that a more vigorous 
entrainment process occurred farther upstream. Towards 
the end of the measurement range, however, the jet cross 
section in this case has become approximately round (Fig. 
3) and the slope of the flux curve has become comparable 
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to that of the corresponding no tab case. Thirds as with 
the subsonic case the rectangular jet with two tabs is 
found to involve the most vigorous entrainment within the 
measurement range. Fourth , the flux values for the lobed 
nozzle are found to be even lower than those for the 
circular nozzle! Finally , one finds that for both the elliptic 
and the rectangular jets the fluxes are significantly higher 
than the circular jet. This is in contrast to the subsonic 
case where the corresponding difference was smaller. 
These points are further discussed next. 

3 . 1 Asymptotic entrainment rate for initially compressible 
jets : 

In regards to the lower entrainment rate for the 
circular jet at the supersonic condition, the following can 
be said based on similarity reasonings. Ideally, all jets 
under consideration, with any nozzle geometry, subsonic 
or supersonic, should eventually become round and self- 
similar with incompressible flow. When that happens, a 
local length scale (Dh) should vary linearly as x and a 
local velocity scale (Uc) inversely as x. 20 It follows that 
the streamwise mass flux m should equal to p A Ci.U c .D H 2 , 
which should vary linearly as x; here p A is the ambient 
fluid density and C t a constant. For the asymptotic 
region, the functional relationships for D H and U c can be 
written as: D H /D = ^.(x/D), and U C /U E = (^/(x/D); 
where C* and C3 are constants. It follows from these 
expressions that the entrainment rate, d(m/m£)/d(x/D), 
equals a constant, C’= Q.C/.Q.^/it). Experimental 
data for incompressible jets yield a value of this constant 
in the approximate range of 0.26 - 0.3. 19,21 ’ 23 

For the asymptotic region of initially compres- 
sible jets, the functional relationships for D H /D and U C /U E 
should still hold. It then follows, 3(m/n%)/3(x/D) = 
(Pa/P^Ci-^ Q-C 4 ^)* Th e constant on the right hand 
side here differs from that for the incompressible case by 
a factor of (p a /Pe)> which for Mj = 1.63 turns out to be 
about 0.65. This qualitatively agrees with the lower slope 
of 0. 18 at the supersonic condition compared to the value 
0.27 observed for the subsonic condition. Here, it is 
pertinent to mention the work of Ref. 24. The functional 
dependence of U C /U E on x/D was investigated based on 
data from earlier experiments. Correlation of the data 
exhibited a mild dependence of U C /U E also on (p A / P e). 
Such a dependence would yield a different expression for 
the entrainment rate. However, the data examined in Ref. 
24 covered the near flow fields. It appears that the far 


asymptotic regions of especially high Mach number jets 
were not covered in the investigation. 

The validity of the expression for the asymptotic 
entrainment rate for initially compressible jets was further 
examined with the axisymmetric nozzle. Normalized mass 
flux curves as in Fig. 7 were obtained for different Mach 
numbers (M,). These data are presented in Fig. 8. The 
slopes calculated from the farthest two data points in each 
curve, and the resultant constants in the expression for 
9(m/ni£)/3(x/D), are listed in Table 1. It is clear that the 
value of the constant, after the ’compressiblity effect’ is 
taken into account, as shown in the third column, is 
within reasonable scatter. Thus, the expression provided 
in the previous paragraph for the asymptotic entrainment 
rate appears to be valid. Note that the rate for Mj = 1.38 
case shows the largest deviation. As to be discussed in 
§3.3, the jet at M x = 1.38 involves screech with a 
pronounced flapping motion. The jet cross section, even 
on the time-averaged field, is somewhat noncircular which 
persists within the measurement range. Thus, as with the 
asymmetric cases (Fig. 7), the jet in this case may not 
have reached the asymptotic state and the higher entrain- 
ment rate is likely to be real. 



Fig. 8 For the circular nozzle, streamwise variations of 
normalized mass flux at different Mj (letters in paren- 
theses indicate screech stages). 
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Table 1 Entrainment rates for the various cases of Fig. 8 


M, 

(Pa/Pe)*C’ 

C’ 

0.30 

0.274 

0.279 

0.65 

0.241 

0.261 

0.95 

0.229 

0.271 

1.14 

0.220 

0.277 

1.38 

0.251 

0.347 

1.63 

0.176 

0.270 

1.80 

0.180 

0.296 

1.96 

0.134 

0.238 


Note, again, that the inferences made in the 
foregoing pertain to the asymptotic regions of the jets. 
Some of the asymmetric jets (Figs. 6 and 7) obviously 
have not reached that state within the measurement range. 
It should also be noted that the effect considered here is 
not the same as the well known coiqpressibility effect 
causing a reduction in mixing layer spreading. 10 The latter 
effect would prolong the length of the developing region 
of the jet. Here, the lower slope measured in the asymp- 
totic region occurs simply due to the fact that the nor- 
malizing parameter, ri^, is relatively large for the initially 
compressible case because of higher ini tial densities. 
Furthermore, it is worth noting in passing that the effort 
to increase jet spreading using various techniques, such as 
tabs, asymmetric nozzle shaping, screech, multiple jet 
resonance, etc., pertain only to the developing regions. 
Far enough downstream, the same asymptotic behavior 
may be expected in all cases. However, it is clear that the 
developing region can be very significantly affected by 
these different techniques which is of great importance in 
practical applications. 

3.2 Effect of tabs and streamwise vorticity on spreading: 

The second and third points made with Fig. 7 
concern the increased entrainment caused by the tabs. The 
observed increase is due to the action of streamwise 
vortex pairs produced by the tabs. This becomes clear 
from the results presented in Refs. 12 and 14, which are 
briefly summarized here. 

A tab produces a pair of counter-rotating stream- 
wise vortices with a sense of rotation such that, between 
the two vortices, ambient fluid is ingested into the core of 
the jet. This is referred to as an ’in-flow’ vortex pair. 
Four equally spaced tabs on the circular nozzle, therefore. 


produce four such in-flow pairs. Vortex induced velocities 
at first make all four pairs move towards the jet axis. 
However, this modem is restricted and, shortly down- 
stream, vortices from adjacent pairs form four ’out-flow’ 
pairs. Each of these pairs causes an ejection of jet core 
fluid into the ambient. The induced motion also propel the 
vortex pairs themselves away from the jet axis. The vigor 
of the dynamics depends on the strength of the vortices 
which in turn depends on the tab geometry. For the 
present tab configuration, the dynamics are indeed quite 
vigorous which cause the observed increase in jet spread- 
ing. In Ref. 12, it was also shown that for the given size 
and geometry, four tabs increased the spreading optimal- 
ly. When six tabs were used, for example, amalgamation 
among the vortices occurred, dropping the number of out- 
flow pairs and reducing the resultant jet spreading. 

With the asymmetric jets the tab effect is more 
complex. This is because these jets typically already 
contain streamwise vortices due to secondary flow within 
the nozzle. Depending on the sense and strength of these, 
and the vortices produced by the tabs, the interaction can 
be different. Let us examine the streamwise vorticity 
distribution for the present tab case with the rectangular 
nozzle. Hot wire measurements were conducted with a 
nozzle which was geometrically similar to the present 
rectangular one but about five times larger in D. The 
larger size was necessary from probe resolution con- 
siderations. Two X-wire probes were used to obtain the 
distributions of the lateral velocities (v and w) on a cross 
sectional plane, the gradients of which provided stream- 
wise vorticity (<*>*). Further details of the measurements 
rm be found in Ref. 14 and those are not repeated here 
for brevity. 

The vorticity data for the tab case are compared 
with the no tab case in Fig. 9. Two iso-surfaces ((oJD/U E 
= ±0.08) are shown. Looking from downstream, darker 
and lighter surfaces represent clockwise and anti-clock- 
wise sense of rotation, respectively. The flowfield without 
tabs can be seen to contain two out-flow vortex pairs 
occurring on the narrow edges. The two pairs would 
move away from the jet axis, and thus tend to keep the jet 
cross section elongated in the major axis plane. By x/D = 
8, the range covered in Fig. 9, these vortices dimmish in 
strength and the jet cross section becomes almost round 
(Fig. 2c). When the tabs are applied, two strong in-flow 
pairs are produced on the narrow edges. These overwhelm 
the naturally occurring vortices and completely dominate 
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fig. 9 Time-averaged streamwise vorticity, o) x D/U E - 
±0.08 iso-surfaces, for a 3:1 rectangular nozzle: (a) no 
tab, (b) two tabs. 


the flow. The in-flow pairs first move towards the jet 
axis, and thus, cause the rapid axis switching. Thereafter, 
the vortices rearrange to form two out-flow pairs which 
continue to cause the lateral spreading in the direction of 
the minor axis. It should be mentioned here that the axis 
switching could be achieved with much smaller size tabs, 
(or prevented by putting the tabs on the long edges; see 
Ref. 14). However, the enormous lateral spreading, as ob- 
served in Figs. 2, 3, occurred only when the tab base 
spanned the entire narrow edge of the nozzle. 

Thus, it is amply clear that the effect of the tabs 
occurs primarily due to the dynamics of streamwise vortex 
pairs. Even though the tabs produce in-flow pairs, the 
vortices rearrange to form out-flow pairs which cause the 
spreading. For the asymmetric jets the process also readily 
affects the axis switching phenomenon. It appears that 
generating an op timum number of out-flow vortex pairs 
with appropriate strength may hold a key to achieving 
even larger jet spreading. 


Note that the dynamics discussed in the foregoing 
pertain to the time-averaged vorticity field. Instantaneous- 
ly the vortices can be much stronger and the dynamics 
more vigorous. 14 The reorientation of azimuthal vorticity 
gives rise to streamwise vortex pairs in the unsteady field 
which also contribute to jet spreading. Liepmann and 
Gharib 25 convincingly demonstrated the important role 
played by instantaneous streamwise vortices in the 
ftntrainm ent process. The role of azimuthal vorticity, on 
the other hand, especially in the near field of an asym- 
metric jet, can also be quite significant for jet spreading 
as well as axis switching, and this is addressed in the next 
section. 

3.3 Effect of screech on jet spreading: 

The fourth point made with Fig. 7 was that the 
least amount of spreading was exhibited by the lobed 
nozzle. An underlying concept for a lobed nozzle is to 
stretch the mixing layer region exposed to the ambient so 
that there is more mixing and entrainment. Thus, the 
result that the lobed nozzle produced less jet spreading 
compared to even the circular nozzle came as a total 
surprise. Soon the reason became apparent. At Mj = 
1.63, all the jets (without tabs) involved screech, whereas 
there was no screech with the lobed nozzle. It may be 
reasoned that the presence of screech with the other 
nozzles and the absence thereof with the lobed nozzle 
caused the difference. 

First, the far field noise spectra for the six cases 
under consideration are shown in Fig. 10. These data 
were obtained with a microphone located approximately 
70D away from the jet axis and about 110° from the 
downstream axis. For ease of comparison the data are 
shown in three pairs. The fundamental screech component 
at about 9 kHz can be clearly seen for the circular jet. 
This is eliminated when the tabs are applied. Both the 
rectangular and the elliptic nozzles are also characterized 
by screech noise. Use of the tabs does not completely 
eliminate screech from the rectangular nozzle but the 
amplitudes are reduced drastically together with a shift in 
the frequencies. It is also clear that the lobed nozzle does 
not involve any screech noise. This is true throughout the 
Mach number range covered in the study. 

Screech is a phenomenon that occurs in imper- 
fectly expanded jets due to the interaction of coherent 
azimuthal vortical structures with the standing shock- 
expansion cells. 26 * 30 The vortical structures are organized 
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% 10 Far-fidd noise spectra for the six cases; the 
three pairs of traces are staggerred successively by one 
major ordinate division. 

and rendered periodic through a feedback loop. The 
complete mechanism remains unclear, but the interaction 
also results in the emission of a tone "...of definite 
frequency and high intensity". 27 

Screech is also accompanied by an increased jet 
spreading. The first known investigation addressing this 
effect for a circular jet is that of Glass. 31 The generation 
of screech in Glass’s study caused a ’precipitous drop’ in 
the ’impact pressure’ measured on the jet axis, clearly 
indicating a significant increase in the spreading of the 
jet. Later on, the effect was further investigated by 
Sherman et al. 32 , and, for a rectangular jet, by Krothapal li 
et al. 33 

Based on past studies cm the effect of artificial 
excitation on incompressible jets, the following comments 
can be made on the likely mechanism for the effect of 
screech on jet spreading. A primary effect of artificial 
excitation, when imparted at an appropriate frequency, is 
to organize and intensify the coherent vortical structures. 
For the case of a circular jet, the excitation causes the 
azimuthal vorticity to periodically roll up into toroidal 
structures. These structures enlarge as they convect 


downstream, various interactions can take place, until they 
break down into disorganized structures near the end of 
the potential core. The process usually is accompanied by 
an increased entrainment. For example, in Ref. 19 a large 
increase in the entrainment was observed under acoustic 
excitation, imparted at the ’preferred mode’ Strouhal 
number (fKcrrATioND/UE = 0.3). About 15% increase in 
the volume flux was measured in the x/D range of 8-10. 
The additional entrainment occurred mostly near the end 
of the potential core (x/D - 4) where the intensified 
vortical structures started to break down. It is the or- 
ganization flud intensifi cation of the coherent structures 
under the screeching condition which are also believed to 
lead to increased jet spreading. The interaction and 
breakdown of the azimuthal vortical structures, organized 
and intensified, in this case, by the self excitation, is 
thought to cause the increased entrainment. 

The relatively higher entrainment with the elliptic 
and rectangular jets at supersonic condition, compared to 
that at subsonic condition, the last observation made with 
Fig. 7, is also believed to be due to screech. Periodic 
excitation can have a more profound effect on an asym- 
metric jet than on an axisymmetric jet. In Ref. 2, acoustic 
excitation of subsonic elliptic jets has been shown to cause 
a dramatic increase in the jet spreading. Under the 
excitation there is an organization and intensification of 
the rolled up asymmetric vortical structures. The self 
induction and the resultant sequence of contortion of these 
structures not only cause axis switching, 2,14 but also result 
in an increased engulfment of ambient fluid. Again, it is 
thought that a similar mechanism is in play with the 
asymmetric jets when screech occurs. This causes the 
larger increase in the fluxes at the supersonic condition 
(Fig. 7) compared to the increase at the subsonic con- 
dition (Fig. 6). This also explains the axis switching seen 
in the supersonic cases (Fig. 3) which did not occur in the 
subsonic cases (Fig. 2). 

Recall the difference in the mass flux results 
between the present asymmetric jets (Fig. 2) and that of 
Ref. 1, as discussed earlier. It should be recognized that 
the dynamics of the asymmetric vortex structures des- 
cribed in the previous paragraph is also active in a natural 
jet (without excitation) because there is natural roll up of 
the structures. Excitation accentuates the process and 
thereby increases the entrainment as well as causes the 
earlier axis switching. The process in the natural case, 
however, can be vary from jet to jet because there can be 
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varying amo unts of self excitation d^ieodisg on the initi al 
and background conditions. At the lower Mach number 
and Reynolds number, die initial boundary layer was 
laminar and the core turbulence was lower in the jets of 
Ref. 1. At M, = 0.3 with the present nozzles, the initial 
boundary layers were nominally turbulent. 14 Thus, a more 
organized initial roll up of the azimuthal vorticity in the 
cited work was likely. This is thought to have ca u se d the 
larger entrainment , in a similar manner as it occurs under 
artificial excitation. 

The discussion so far has been based on previous 
involving ’plane wave’ acoustic excitation. There 
have been other studies involving excitation with more 
complex waveforms (subharmonic, helical, combination 
of modes, etc.). In most cases, artificial excitation at the 
various modes also cause an increase in the overall jet 
spreading. 34 Screech involves self excitation in a variety 
of modes. The characteristics of these modes and their 
impact on jet spreading are further explored in the fol- 
lowing. 

3.3.1 Unsteady flow characteristics in various stages of 
screech: First, let us examine the screech frequency 
chracteristics. Variations of screech frequency, expressed 
as Strouhal number, for the circular, elliptic and rectan- 
gular jets are shown in Fig. 11. The triangular and 



Kg. 11 Nondimensional screech frequency (fj) versus 
Mj for the indicated nozzles; letters indicate screech 
stages. 


ftlli priral symbols represent the rectangular and elliptic 
jets, respectively, while the dotted line represent data for 
the circular jet. As observed by many previous resear- 
chers, 27 ' 29 ’ 35 ’ 36 the screech for the circular nozzle is seen to 
go through different stages as the nozzle pressure ratio (or 
Mj) is varied. These stages involve, as noted in the cited 
studies ar| d as to be further discussed in the following, 
flow field o s cillations in the axisymmetric (stages A, and 
Aj), flap pin g (B and D) and helical (C) modes. Corre- 
sponding frequency characteristics for the rectangular and 
the elli ptic jets also show a staging behavior involving 
one frequency jump within the measurement range. A 
jump occurs from B’ mode to D’ mode around Mj = 1.65 
(52 psig) for both the nozzles, with a certain amount of 
overlap between modes. The nature of the unsteady flow 
fields in the different screech stages for the different 
nozzles is further explored by the phase averaged total 
pressure measurements. 

As described in §2, the phase averaged data were 
obtained at a given downstream location over the jet cross 
sectional plane. The location of measurement was deter- 
mined for each Mj for a given nozzle by initial surveys. 
The choice was dictated by the requirements that the 
location was far enough downstream in order to avoid 
probe interference as much as possible, yet sufficiently 
upstream to ensure coherence in the flow fluctuations. 
These criteria yielded measurement locations that were 
generally farther downstream with increasing M,. Typical- 
ly, the data were acquired for 13x13 grid points over the 
cross sectional plane. The signal from a microphone 
located near the nozzle exit was used as reference. The 
phase averaged total pressure, relative to the time ave- 
raged total pressure, was measured at each grid point. The 
were obtained for nineteen phase points over a cycle 
and involved averaging over approximately 300 sample 
functions. 

Data for stages A l9 A 2 , B, C and D of the 
circular nozzle are shown in Fig. 12. The Mach numbers 
and measurement locations are indicated in the figure 
caption. For each case, two equal and opposite (positive 
and negative) iso-surfaces are shown. The viewing angle 
(same for all the cases) and the iso-surface levels were 
chosen to demonstrate the nature of the fluctuating field 
as clearly as possible. The data vividly demonstrate that 
the mode shapes in the two A stages are axisymmetric. 
Those in stages C and D are clearly helical and flapping, 
respectively. The flow pattern in stage B is also primarily 
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Fig. 12 Phase-averaged unsteady total pressure dis- 
tributions for the circular jet. Screech stage, Mj and 
measurement location (x/D) are: (a) A lf 1.11, 4.5, 

(b) A 2 , 1.16, 4.5, (c) B, 1.37, 6, (d) C, 1.63, 9, 

(e) D, 1.81, 12. 


flapping but is apparently mixed with a helical pattern. 
An unsteadiness and switching from helical to flapping 
oscillation was also inferred for the B mode in the studies 
of Refs. 35 and 37. As stated before the mode shapes 
associated with the circular jet screech stages had been 
inferred earlier in several studies. These were done 
through correlation measurements in the near acoustic 
field as well as by schlieren visualization of the flowfield. 
However, ambiguities can arise in interpreting the results 
of such measurements. For example, it can be difficult to 
differentiate the helical mode from the flapping mode. 
The present results unambiguously confirm the mode 
shapes and provide a visual impression of the fluctuating 
flowfield. 

The corresponding data for the stages denoted B’ 
and D’ (Fig. 11), for the 3:1 rectangular nozzle, are 
shown in Fig. 13. It is clear that both the stages involve 
flapping motion about the x,y plane. Similar data for the 
elliptic nozzle (not shown) also yielded essentially the 
same modal structures. 



Why does the jet go through a frequency jump 
from stages B’ to D’ with seemingly similar unsteady 
motion? Why does the circular jet go through a jump 
from Aj to A 2 mode both having axisymmetric unsteady 
motion? What dictates the other modal structures with the 
circular nozzle? What chooses the direction of the helical 
motion as well as the plane about which the flapping 
motion occurs? These issues, as well as the complete 
mechanism for the screech noise generation, remain far 
from clear. These are considered beyond the scope of the 
present work. (The reader is referred to Ref. 36 for an 
insight into the unsteady shock motion associated with 
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Fig. 13 Phase-averaged unsteady total pressure dis- 
tributions for the rectangular jet. Screech stage, Mj 
and x/D are: (a) B’, 1.63, 6, (b) D’, 1.75, 8. 

screech, and to Ref. 38 for the effect of induced screech 
on the jet via placement of obstacles in the flow. These 
are parallel studies being conducted in the same labora- 
tory). Here, let us make some further comments on the 
data of Figs. 12 and 13 before turning attention back to 
the spreading of the various jets. 

The unsteady amplitudes shown in Figs. 12 and 
13 are values over the local time averaged values. The 
maximum peak-to-peak amplitudes are typically about 
0.2% of the local total pressure. These are small com- 
pared to common experience with excitation studies of 
subsonic jets. The red (darker) iso-surface represents 
regions of positive (higher) total pressure or streamwise 
velocity. The yellow (lighter) iso-surface represents 
regions of negative (lower) streamwise velocity. For the 
axisymmetric modes (A t , with the circular jet, one 
would expect the toroidal vortical structure to be located 
near the middle of the field, where the velocities are 


higher in the core but lower outside. (Note that an inner 
red iso-surface in this region is surrounded by a outer 
yellow iso-surface). Mode C (Fig. 12d) is clearly helical. 
The location of the corresponding helical vortex tube is 
not completely clear, but is probably between the red and 
yellow regions. Mode D can be viewed as due to the 
superposition of two opposite sense helices causing the 
flapping motion. A certain plane of flapping is chosen by 
the jet, presumably, either due to irregularities in the 
geometry of the upstream hardware or in the recirculating 
flow in the laboratory. For the rectangular nozzle, the 
flapping plane is well defined and coincident with the 
major axis plane of the nozzle. This is clearly seen in 
both sets of data of Fig. 13. 

3.3.2 Effect of various screech stages on jet spreading : 
The normalized mass flux variations, measured at x/D = 
14, as a function of Mj for the circular jet are shown in 
Fig. 14. The corresponding screech stages (Fig. 1 1) are 
reproduced on the top for comparison. One finds that the 
fluxes undergo large variations with Mj. The fluxes are 
highest in the flapping mode B, and are again relatively 
high in the next flapping mode D. With the onset of the 
helical mode C, the flux values drop substantially. Pitot 
probe data, obtained at a single point on the jet centerline, 
reported in Refs. 31 and 32 are in general agreement with 
these results. For example, the total pressure was ob- 
served in those studies to be lower at pressure ratios 
corresponding to the B mode indicating an increased jet 
spreading. 

Referring back to Fig. 7, it is reasonable to infer 
that without screech the fluxes for the circular jet would 
have been lower than those for the lobed nozzle. The 
occurrence of screech generally results in an increased jet 
spreading and this effect is most pronounced in the 
flapping modes. Note that the data in Fig. 7 for the 
circular jet (at Mj = 1.63) involved screech in the C 
mode which is characterized by relatively less spreading. 
Therefore, the difference between the circular and the 
lobed nozzle data would have been even more dramatic at 
either a lower or a higher M, yielding the B or D screech 
modes for the circular nozzle. 

Corresponding mass flux data, measured at x/D 
= 14, for all six cases are compared in Fig. 15. The 
fluxes for the elliptic and the rectangular jets are clearly 
high in the B’ stage. The pattern is similar to the circular 
jet B stage. However, it is not clear why the flux values 
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fig, 14 Normalized mass flux measured at x/D = 14, 
as a function of Mj, for the circular nozzle* Nondimen- 
sional screech frequency (Fig. 11) is shown on the top* 


drop as the D’ mode, which is also a flapping mode, is 
approached. It is possible that the screech amplitude, 
which diminishes as a stage jump is approached, is 
responsible for the trend. In general, however, it is clear 
that the flapping mode screech results in the highest jet 
spreading. This is consistent with the acoustic excitation 
results of Ref. 39 in which a combination of + 1 and -1 
modes (opposite sign helices) also yielded a remarkable jet 
spreading. 

Referring back to the discussion earlier in this 
section, an additional reason for the relatively higher 
fluxes for the asymmetric jets, seen in Fig. 7, becomes 
clear. At Mj = 1.63, represented in Fig. 7, the elliptic 
and the rectangular jets wait through a flapping mode 
screech whereas the circular jet wait through a helical 
mode screech. This contributed to the observed higher 
fluxes for the asymmetric jets. 

Once again, let us note that the fluxes for the tab 
cases in Fig. 15 are much larger throughout the Mach 
number range. The values are larger than those obtained 
with any of the asymmetric nozzles. This is true for 
supersonic as well as subsonic conditions. Note that 
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Kg. is Normalized mass flux measured at x/D — 14, 
as a function of Mj, for all six cases. 


especially the circular jet with tabs does not involve 
screech. Whereas the spreading increase with screech 
occurs through an organization of the azimuthal vortical 
structures and their subsequent dynamics, that with the 
tabs occurs primarily through the dynamics of the stream- 
wise vortex pairs. 

3.4 Thrust loss for the tab cases: 

The measured thrust as a function of the nozzle 
pressure ratio is shown in Fig. 16 for the six cases. The 
data were obtained by a simple device involving a Toad 
cell 9 and an auxiliary plenum mounted on linear bearings, 
details of which can be found in Ref. 12. The data are 
shown in two sets, the lower set is for the circular jet and 
the upper set is for the other cases. The upper set is 
shifted by one major ordinate division. Since the circular 
jet had a smaller diameter, the corresponding thrust values 
are Iowa. In each set the ideal thrust is shown by the 
solid line. It becomes apparent that all the no tab cases 
involve slightly lower thrusts relative to the ideal value. 
The somewhat Iowa values are due to the boundary 
layers within the nozzles. The ideal thrust is calculated 
with the assumption of uniform flow with zero boundary 
layer thickness. 
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When the tabs are added, as expected, there is 
considerable thrust loss. For example, at Mj = 1.63 
(P T /p A = 4.47), the thrust loss for either of the two tab 
cases in Fig. 16 was approximately 12%. However, with 
the tabs in place, there is area blockage and the mass flow 
is lower for a given P T /p A . An appropriate evaluation of 
the performance of a given nozzle configuration is 
obtained by measurement of the thrust coefficient, which, 
for the six cases, is shown in Fig. 17. These data required 
simultaneous measurement of the actual mass flow which 
was done by an orifice meter located in the supply line. 
Unfortunately, the uncertainty in the mass flow data was 
relatively large and resulted in a scatter of about ±2% in 
the thrust coefficient. 

It should also be mentioned here that the ideal 
thrust, used to obtain C f , was calculated by assuming 
convergent flow (which was the case). For a given P T /p A , 
full expansion (via convergent-divergent nozzle) would 
yield the maximum available thrust. 40 This would be 
higher, for example, by about 2% at P T /p A = 4.47. If this 
were used in the calculation, the thrust coefficients would 
have correspondingly lower values. Note also that the data 



fig. 17 Thrust coefficient versus nozzle pressure ratio 
for the six cases. 

in Fig. 17 are shown only in a high range of P T /p A where 
the measurement uncertainty was lower and within the 
limits stated in the foregoing. Least-squares-fitted curves 
through approximately six data points are shown. 

As expected, the thrust coefficient turns out to be 
nearly unity for all the no-tab cases throughout the nozzle 
pressure ratio range covered in Fig. 17. With the tabs, of 
course, there is a performance penalty. However, the 
penalty is only about 5 % for the rectangular jet with two 
tabs, and about 7% for the circular jet with four tabs. 
Recall that the corresponding geometric area blockages 
were about 12% and 8%, respectively. Thus, the per- 
centage penalty for the circular jet case turns out to be 
about the same as the area blockage. The corresponding 
penalty for the rectangular case, on the other hand, is 
lower. The difference borders on the uncertainty in the 
data but appears to be real. Unfortunately, not much can 
be said about the fluid dynamics from the thrust data, and 
an explanation for the difference will require further 
investigation of the flow field close to the tabs. It should 
be pointed out that in the case of four tabs, as opposed to 
two, the flow has to negotiate more comers. Also, with 
the circular nozzle, fabrication resulted in a small part of 
the tab, (sector determined by the base of the tab and the 
curved nozzle wall), standing perpendicular to the flow. 
In the rectangular nozzle case, the entire tab surface is at 
45° to the stream direction. These factors could explain 
the relatively higher performance penalty in the tab case 
with the circular nozzle. 
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Concluding Remarks: 

The spreading characteristics of free jets from 
various asymmetric nozzles have been studied. The 
nozzles included a circular, a 3:1 rectangular, a 3:1 
elliptic and a six-lobed one. In addition, jets from the 
rectangular nozzle fitted with two tabs on the narrow 
edges and the circular nozzle fitted with four equally 
spaced tabs were also included in the comparative study. 
Measurements covered incompressible subsonic conditions 
to underexpanded supersonic conditions up to a jet Mach 
number of 1.96. 

Measurements in the asymptotic regions of the 
circular jet show that the entrainment rate, when nor- 
malized by the initial mass flux, the nozzle diameter and 
the ratio of the density between the ambient and the 
nozzle exit, is a constant. The constant has a value of 
approximately 0.27, agreeing with previous results for 
incompressible jets. Most asymmetric jets, however, did 
not reach the asymptotic state within the measurement 
range covering thirty equivalent diameters from the nozzle 
exit. Screech, nozzle shape and the tabs altered the jet 
spreading within the developing regions very signifi- 
cantly. 

At subsonic conditions the spreading of die 
asymmetric jets was found to be only slightly more than 
that of the circular jet. But it was found to be signifi- 
cantly more at supersonic conditions when screech 
occurred. Screech increased the jet spreading with all 
nozzles, the amount of increase varying with varying 
stages of the screech. The lobed nozzle which did not 
involve screech, therefore, exhibited the least spreading at 
supersonic conditions. 

The unsteady flow characteristics under the 
various stages of screech, for the circular and the rectan- 
gular nozzle have been studied via phase averaged measur- 
ement of the fluctuating total pressure. The axisymmetric, 
flapping and helical unsteady motions occurring in the 
vaious stages of screech with the circular nozzle are 
vividly demonstrated by these data. Two stages of screech 
are observed with the rectangular and the elliptic nozzles. 
Both stages involve the flapping unsteady motion. Among 
the various stages, screech involving the flapping motion 
caused the maximum jet spreading. 

Jet spreading for the two tab configurations was 
found to be the largest among all the cases. This was true 
in the subsonic regime, as well as in the supersonic 
regime in spite of the fact that screech was eliminated by 


the tabs. The dynamics of the streamwise vortex pairs 
produced by the tabs caused the most efficient jet spread- 
ing. Thus, a manipulation of the latter dynamics is likely 
to hold the key to possible further increases in jet spread- 
ing needed in various applications. 
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